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Abstract 

Schwinger-Dyson equations for propagators are solved for the scalar 
theory and massive Wick-Cutkosky model. With the help of integral represen¬ 
tation the results are obtained directly in Minkowski space in and beyond bare 
vertex approximation. Various renormalization scheme are employed which 
differ by the finite strength field renormalization function Z . The S'—matrix 
is puzzled from the Green’s function and the effect of truncation of the DSEs 
is studied. Independently on the approximation the numerical solution breaks 
down for certain critical value of the coupling constant, for which the on-shell 
renormalized propagator starts to develop the unphysical singularity at very 
high space-like square of momenta. 
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1 Introduction 


The Dyson-Schwinger equations (DSE) are an infinite tower of coupled integral equa¬ 
tion relating Green functions of the quantum field theory. If solved exactly, they 
would provide solutions of the underlying quantum field theory. In practice, the sys¬ 
tem of equations is truncated and one hopes to get some information, in particular on 
the solution in the non-perturbative regime, from solving the simplest equation for 
the two-point Green functions- the propagators. The other vertex functions, which 
also enter the DSE for the propagator, are either taken in their bare form or some 
physically motivated Ansatze is employed. 

In most papers dealing with the solution of DSEs, the Wick rotation from the 
Minkowski to Euclidean space is used to avoid singularities of the kernel inherent to 
the physical Greens functions. To our knowledge, the only exception is the series 
of papers [1],[2],[3],[4] employing so-called ’’gauge technique” in quantum electrody¬ 
namic and its gauge invariant extension to quantum chromodynamic [5], this work 
represents the born of the ”Pinch Technique”. Until now, the above-mentioned ap¬ 
proach has been never used in its non-perturbative context. Although not dealing 
with gauge theories, similarly to these techniques we instead solve the directly in the 
momentum space, making use of the known analytical structure of the propagator, 
expressed via the spectral decomposition. In the spectral or dispersive technique, we 
write the Green function as spectral integral over certain weight function and denom¬ 
inator parameterizing known or assumed analytical structure. The generic spectral 
decomposition of the renormalized propagator reads: 

G(p 2 ) = f da a(a) . , (1) 

J p 2 — a — le 

where a(a) is called Lehmann weight or simply the spectral function. If the threshold 
is situated above the particle mass, as it is for the stable (and unconfined) particles 
then the spectral function typically looks like 

a(a) = rS(m 2 — a) + a (a) (2) 

where the singular delta function corresponds with non-interacting fields and cr is 
appeared due to the interaction. Finite parameters r then represents the propagator 
residuum and is simply related to field renormalization. It is also supposed that a is 
a positive regular function which is spread smoothly from the zero at the threshold. 
Note here that the positivity of Lehmann weight is not required for our solution, but 
the models studied in this paper naturally embodied this property ;see for instance 
[6], or any standard textbook. 

Putting the spectral decomposition of the propagators and the expression for 
the vertex function into the DSE allows one to derive the real integral equation 
for the weight function cr(a). This equation involves only one real principal value 
integration and can be solved numerically by iterations. Our solutions are obtained 
both for space-like and time-like propagator momenta; obtaining this in the Euclidean 
approach would require tricky backward analytical continuation. Since all momentum 
integration are performed analytically, there is no numerical uncertainty following 
from the renormalization which is usually present in Euclidean formalism [7]. Here, 
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the renormalization procedure is performed analytically with the help of the direct 
subtraction in momentum space. This perturbative perfectly known renormalization 
scheme (see [8] for scalar models and [9] for QED case, where a the comparison with 
other perturbative renormalization schemes was also made) has been already applied 
to the QED and Yukawa model [10] in its non-perturbative context. In this paper, the 
off-shell momentum subtraction renormalization scheme was introduced and used. In 
order to simplify the technique and to compare various schemes we restrict ourselves 
to the choice of on mass-shell subtraction point /i = m. 

In this work, we would like to present certain solutions of rather obscure theories: 
<f> 3 and 'hf'hj scalar models. The second model will be referred here as the (gen¬ 
eralized) Wick-Cutkosky model (WCM). In fact, not only models mentioned above 
but the all super-renormalizable four-dimensional scalar models are not properly de¬ 
fined since they have no true vacuum [11], Instead of this they have only metastable 
vacua (here we assume non-zero masses of all particle content, in the opposite case 
the appropriate classical potential would not posses any local minimum). Instead of 
discarding these types of models, as sometimes happens, we look whether this ’in¬ 
consistency’ can be captured by the formalism of DSEs, or whether the appropriate 
solutions ’behaves ordinarily’. The property of super-renormalizability makes our 
models particularly suitable for this purpose. Actually, the super-renormalizability 
here implies the finiteness of the renormalization held constant Z which therefore can 
not be considered at all (i.e.Z = 1). In the case of <h 3 theory we do not fully omit 
the held renormalization ,but with the help of the appropriate choice of the constant 
Z, we choose the given renormalization scheme. Making this explicitly and after 
the evaluation of the scattering amplitude we look (in each scheme) at whether the 
observables do converge (in all schemes) to some experimentally measurable values 
of the virtual scalar world . As the suitable observables we choose the amplitude M 
for the scattering process $$ —» $$ and we have no find any unexpected or even 
pathology behavior. Instead of this, when the approximation of the full solutions 
improve, we will see that the amplitudes calculated in the various renormalization 
schemes tends to converge to each other, i.e. in this aspect, the <h 3 theory behaves 
as the ordinary and physically meaningful one. Here, this is the right place to note 
that the models with the metastable ground state serve as an useful methodologi¬ 
cal tool, the role in which they are often employed. In fact, <h 3 theory serves as a 
good ground for the study of the various phenomena [12], [13], [14],[15] (including 
phenomena like non-perturbative asymptotic freedom and non-perturbative renor¬ 
malization). There also exist a number of papers dealing with WCM. The DSEs 
for propagators of the WCM in their simple bare vertex approximation have been 
solved for the purpose of calculation relativistic bound states [16] (for other recent 
work dealing with the bound states problem within the WCM see [17] and references 
therein). For the purpose of comparison with [16] we solve the exactly analogical 
Minkowski problem. The obtained value of the critical coupling should depends on 
the renormalization scheme. Having this slight dependence under control it allows 
us to compare with other non-perturbative method [18],[19]. The comparison with 
conventional perturbation theory is also made. 

Regardless of the facts mentioned above, we are far from concluding that $ 3 model 
is a fully physically satisfactory one, since we do not know anything about the full 
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solution. At this point, the study presented in this paper and the studies of the <f> 3 
model in five [12] and six [13] dimensions are conclusive in a similarly cautious way. 
Probably, a more sophisticated conclusion could be obtained by some lattice study, 
which has not yet done for this purpose. 

At the end of the introduction, we should mention that there is always the pos¬ 
sibility of including a sophisticated cut-off function /(A) into the Lagrangian and 
regard onr cubic models as an effective model bellow this cut-off. The theory at en¬ 
ergies above A could be another field theory or string theory, or whatever. However, 
this method is developed and the appropriate Polchinski renormgroup equations may 
be written down [20], [21] these cut-off methods he somehow beyond the scope of 
this paper and we prefer to use the usual renormalization schemes, where the inde¬ 
pendence on the appropriate regularization procedure is manifest. Clearly, with the 
use of cut-off method it would be difficult to perform the aforementioned comparison 
with the results [16], where on mass shell renormalization scheme has been performed. 
Furthermore, we should note here that the dispersion technique used thorough the 
proposed paper would become more complicated due to the presence of the profile 
function /(A). 

In the next section we present the DSEs for <f> 3 for propagator and vertex func¬ 
tions. Subsequently we discuss the renormalization procedure and rewrite the prop¬ 
agator equation into its spectral form. Also, the numerical results and limitations 
are discussed. The WCM is dealt with in the section 3. It is solved numerically in 
its pure bare vertex approximation. The details of calculations are relegated to the 
appendices A and B. 

2 <f> 3 theory 

2.1 Dyson-Schwinger equation for theory 

The Lagrangian density for this model reads 

£ = ^d fl (p 0 (x)d fl Mx) ~ \rnl(t>l{x) - g 0 <l>l(x) (3) 

where index 0 indicates the unrenormalized quantities. With the help of the 
functional differentiation of the generating functional (for this procedure, see for 
instance [22]) with classical action determined by (3) one gets the following DSE 
(after transforming into the momentum space) for the inverse propagator 

Gq\p 2 ) = P 2 ~ m o ~ n o(p 2 ) 

n 0 (p 2 ) = i3g 0 J ^-jjT 0 (p ~ q,q)Go(p - q)G 0 (q) (4) 

where To is the full irreducible three-point vertex function which satisfies its own DSE 
(10). The integral of Ho is divergent and requires the mass renormalization. Making 
the on-mass-shell subtraction we define renormalized selfenergy n#i 

n ra (p 2 ) = n 0 (p 2 ) -n 0 (m 2 ) (5) 
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where m is the pole ’’physical” mass, given by the equation G 1 (m 2 ) = 0. Defining 
the mass counter-term 

m 2 = mg — 5 m 2 , 6m 2 — II 0 (m 2 ) (6) 

and introducing additional finite renormalization constant 

0o = '/Zfa go = g~ r \ (T) 

Z 2 

we obtain the inverse of the full propagator in term of physical mass 


G 1 (p 2 ) = Z (p 2 — rn 2 ) — Hi(p 2 ) 

Mp 2 ) = z (n 0 (p 2 ) - n 0 (m 2 )) 

ZU 0 (p 2 ) = i3g 2 j 7 ^ 3 i r (p- q,q)G(p- q)G{q) (8) 

where g is a renormalized coupling and the constant Z g corresponds with the renor¬ 
malization of the vertex function, and G represents the renormalized propagator with 
respect to the field strength renormalization, i.e. 

r = Z g T 0 , Go(p 2 ) = ZG{p 2 ) (9) 

We closed the system of our DSEs already at the level of equation for proper 
vertex. Instead of solving the full renormalized DSE for the vertex 

gT(p,l) = 6g + 6gi J -^-T(p,q)G(q)G(l - q)M(q,l,p) (10) 

we approximate the vertex by the first two terms of the appropriate skeleton expansion 

g r(p,l) = 6g + i(6g) 3 J ^-G(q)G(p - q)G(l - q). (11) 

i.e., we approximate the vertex inside the loop by its bare value and the scattering 
matrix M in Eq. (10) is taken in its dressed tree approximation,i.e., M = G. In 
the following, we will call the dressed vertex (DV) or improved approximation for the 
solution of the propagator when the equation (11) is used for obtaining the triplet 
scalar vertex Y and in the same spirit we use the name bare vertex (BV) for such a 
solution where only the bare vertex was used. The improvement of the approximation 
is achieved by the skeleton expansion of the proper Green function where the series 
of DSEs is thrown away. Here, this is done at the level of the triplet vertex. In [23], 
we can see how the problem is becoming more complicated when M is nontrivially 
taken into account. 

The equation for the propagator is solved in BV and DV approximations at each 
renormalization scheme separately. We define these in the following section. 
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2.2 Choosing the scheme 

We assume (or rather we neglect it) that there interaction does not create the bound 
states contributing to the weight function a. We use the name minimal momentum 
subtraction renormalization scheme (MMS) the one where the only mass subtraction 
is used and where the held leaves unrenormalized, i.e. Z — 1. Therefore, we can 
write the spectral decomposition for the propagator and for the self-energy II in the 
following form: 


G(P 2 ) 


Ih (p 2 ) 


p 2 — m 2 


+ 


r d a 2 ° (a) 

Um 2 p z — a — ie 
-1 


\P 


da 


m 2 — IIi(p 2 )j 

p(a)(p 2 — m 2 


4m 2 (a — m 2 ) (p 2 — a + ie) 


( 12 ) 


where np(s) represents the self-energy absorptive part and the threshold value of 
momentum P 2 = 4m 2 is explicitly written. Obviously, in this MMS the propagator 
does not have the pole residue r equal to unity 


}im(p 2 — m 2 )G(p 2 ) = 


d 


J- 2 nAp% 2 =m 2 = 


da 


dp 2 

-p(a 


ril(P )| p 2 =m 2 


= r 


4 m 2 (a — m 2 ) 2 


(13) 


After a simple algebra and taking the imaginary part of equation (12) we arrive at 
the relation between the spectral functions a and p 


a(u) 


rp{oj) 1 a(u)p(a)^^ + a(a)p(u) 

7--9-' • / Ca - 

(cu — m z ) z uj — m z Jim 2 u) — a 


(14) 


where P. denotes principal value integration. 

This is the first from two necessarily coupled equations which we actually solve 
for a given theory. We discuss it in some details since its form depends only on the 
adopted renormalization procedure, not on the actual form of the interaction, nor on 
the approximation employed for the vertex function T in the DSE for the propagator. 
The second equation connecting o and p does depend on the form of vertex. Its 
derivation is more complicated and we deal with it in the Appendix A. 

In some cases, the form (14) is not the most convenient one; for instance, when 
we want to look the bound state spectrum influence causes just by the self-energy 
effect ([24]) . Note the presence of the constant r in the first term on the right-hand 
side, it has to be determined from the relations (13) after each iteration. To get rid 
of this, we define the usual on-shell renormalization scheme with unit residuum (OSR 
scheme) by 


Z = 1 + 5Z] 5Z 


d 

dp 2 


IIl(p ) |p2_, m 2 


5 


(15) 



which gives the standard receipt how to calculate the OSR propagator 


G osr (p 2 ) = 

[p 2 - m 2 - n 2 (p 2 )} 


n 2 (p 2 ) = 

n(p 2 ) - n(m 2 ) - -^n(p 2 )| P 2 =m2 (p 2 - m 2 ) 


n (p 2 ) = 

*3^2 / ( 27r ) 4 r (p q,q)G(p g)G(q) 

(16) 

and subsequently implies the spectral decomposition for G OSR and II 2 


Gosr{p 2 , 

. 1 f°° , 02 (a) 

) — 2 + / da 

p z — rrr Jim 2 p z — a — ic 


n 2 (p 2 ; 

j p 2 (a)(p 2 - m 2 ) 2 

} “ Jim 2 (a- m 2 ) 2 (p 2 -a + ie)' 

(17) 

The relation between a 2 

and p 2 is now derived in the same way as 

before and it 


reads 


o~2 M 


P 2 M 
(ou — m 2 ) 2 




0- 2 (u;)p 2 (a) 
da - 


+ C r 2(«)P2(<^) 

ca — a 


(18) 


Note, that the Eqs. (14),(18) are inequivalent clue to the scheme difference, the 
appropriate dependence of the weights p and p 2 on the coupling constant g and 
g 2 is explicitly written in appendices A and B, respectively. ( Two inequivalent 
renormalization schemes should give the different Green function, but should give 
the same S-matrix). 

At the end of this section, we very briefly discuss dimensional renormalization pre¬ 
scription [25], showing here that it is fully equivalent to MMS to all orders (note that 
the perturbation theory is naturally generated by the coupling constant expansion 
of the DSEs solution). For this purpose we choose the modified minimal subtraction 
MS scheme, noting that any other sort of schemes based on the dimensional regular¬ 
ization method would be treated in the same way. Since the only infinite contribution 
are affected when this renormalization is applied, therefore the contribution with the 
dressed vertex (master diagram and so that) satisfies the unsubtracted dispersion re¬ 
lation while for instance the one loop skeleton self-energy diagram (in a fact the only 
one irreducible contribution which is infinite in four dimension) looks (for space-like 
momenta) like 



18 g 2 
(47r) 2 



m 2 — p 2 x( 1 — x ) 

dt'Hooft 


+ n (p 2 )/ inite 


(19) 


where II finite represents the omitted finite terms which are not affected by dimensional 
renormalization at all (since they are finite to the all orders). 

The inverse of the full propagator reads in this scheme 
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G Ms(P 2 ) =P 2 - - n^ 5 (p 2 ) - n (p 2 )finite- (20) 

Identifying the pole mass by equality Gjf S (p 2 = m 2 ) = 0 we simply arrive to the 
result 

G-J s (p 2 )=p 2 -m 2 p -U 1 (p 2 ) (21) 

where 

m 2 p = m 2 (n t 'Hooft) + n [ lj s (m 2 p ) + n(m 2 )/ inite 
nib 2 ) = nj] s (p 2 ) - ii(p 2 ) fi n it e - ng 5 (mj) - n(m 2 ) /inite ( 22 ) 

Since the pole mass is renormgroup invariant quantity, we see that MS scheme 
exactly corresponds with the one subtraction renormalization scheme,i.e. the MMS. 
Note here, that in renormalizable models such identification is not so straightfor¬ 
ward but always possible [9]. Of course, the appropriate identification is then rather 
complicated. To conclude this section, we can see that the popular renormalization 
prescription like MS or AdS schemes can be ordinarily used in the non-perturbative 
context. At this point we disagree with the opposite statement of the paper [28]. 

2.3 Test of scheme (in-)dependence 

The physical observables should be invariant not only with respect to the choice of 
renormalization scale, but also with respect to the choice of renormalization scheme. 
The first invariance is more then manifest in our approach, since all the quantities 
used here are the renormgroup invariants. The second mentioned invariance is less 
obvious and, in fact, it is clear only for some very simple cases. (The most simple 
case is the tree-level amplitude evaluation, where the residua of the propagators may 
be exactly absorbed into the redefinitions of the coupling constants; but of course, 
in this case the renormalization is not required ). In any reasonable renormalizable 
quantum held theory it is strongly believed that the obtained exact Green functions 
must build the same S'-matrix. In perturbation theory, we usually have several first 
terms of perturbation expansion and we hope that they offer satisfactory description 
of the nature when the ”right” choice of renormalization scheme is made [9],[26]. 
Furthermore, we should be aware that the possible sum of infinite many terms of 
perturbation series should be regarded as an asymptotic one. In fact, the application 
of some sophisticated resummation technique is necessary in that case [22],[27]. 

In DSE treatment we can talk about the level of DSEs system truncation instead 
of a given coupling order. In the text bellow we describe a simple possible procedure 
how to see the improvements of physical observable when it is calculated within the 
improved truncation of DSEs. For this purpose the BV and the DV solutions of 
DSEs in the both MMS and OSR schema are used to compose the same physically 
measurable quantity. 

For our explanation we have explicitly choosen the matrix element M of the clastic 
scattering process —> (jxj) which can be written 
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( 23 ) 


M(s,t, u)= J2 rG(o)r + ... = ]T (6g) 2 G(a) + ... 

a=s,t,u a=s,t,u 

The dots denote the neglected boxes and crossed boxes contributions and the letters 
s,t,u in (23) represent the usual Mandelstam variables that satisfies s + t + u = 4m 2 , 
since now, the external particles are on-shell. 

Using the notations introduced in the previous section, then the matrix M in 
MMS scheme is calculated as 


M mms = ]T (6 g) 2 G{a), (24) 

a=s,t,u 

where the propagator is calculated through equations (12),(13),(14). For OSR, the 
scattering matrix is composed like 

M OSR = ^2 (6 g2) 2 GosR{a) (25) 

a=s,t,u 

where the propagator is calculated through equations (16),(17) and (18) and the 
relations for p’s are reviewed in the appendices A and B. 

In the ideal case we would obtain 

M M ms = M 0 sr , (26) 

which should be consequence of exact scheme independence. In reality, the Rcl.(26) 
is not exactly fulfilled due to the truncation of DSEs system. In what follows, we 
describe how to check the consistent condition (26) and how to see the appropriate 
deviation numerically. 

Clearly, the equality should be valid in each kinematic channel separately. For 
instance, choosing the f-channel for this purpose and comparing the pole part of 
matrices M we obtain the relation between the coupling at each scheme 

& = rg 2 , (27) 

where r is the residuum calculated from equation (13). This implies for us that if we 
calculate the Green’s functions in the OSR scheme to compose the same S-matrix the 
Green’s functions in the MMS scheme must be calculated with the coupling g 2 = y/rg. 
Having the results for a and <j 2 extracted from the DSEs solved in the appropriate 
schemes, we can compare imaginary parts of scattering matrices Mmms and Mosr- 
Our approximation (24),(25) implies 

= g 2 °(u)- (28) 

How accurately this equality is fulfilled at non-trivial regime t > 4m 2 can be simply 
checked. For this purpose we evaluate the integral (weighted) deviation E N 

g _ f [Mossit) ~ M M Ms(t)\ jw 
I [Mossif] + M MM s(t)\ jw 
S[cr 2 {t)-cr(t)/r}^ 

S[? 2 (t) +a(t)/r] $ 
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( 29 ) 



where the parameter N serves us for adjusting the regime of momenta we are inter¬ 
ested in. A larger value of N enhanced the threshold values of momenta while the 
ultraviolet modes are suppressed in that case. We choose N — 0,1 for the purpose of 
this paper. 

Let us stress at the end of this section, that the next leading order of M is scheme 
invariant and all the difference therefore follow from the remnant of the full DSEs 
solution. Hence only negligible deviation is expected for small couplings. Also in 
general, the deviation E^ should decrease when considered approximations become 
more and more close to the full nonperturbative solution and it should principally 
vanishe for the exact solution. In other words, En must decrease when approximation 
(truncation of DSEs) improves. The results obtained by the above sketched method 
are reviewed in the next section. 

2.4 Results 

The integral equation for Lehmann weights have been solved numerically by the 
method of iteration. The appropriate solutions, obtained for several hundreds of 
mesh points and with the use of some sophisticated integrator, have an accuracy 
of approximately one part of 10 4 for reasonable value (A << A cr u) of the coupling 
strength A and increase (up to several %) when A ~ \ C rit- The coupling strength is 
defined as dimensionless quantity 


1 Sg 2 

167r 2 m 2 


(30) 


The critical value of A is simply defined by the collapse of ( numerically sophisti¬ 
cated) solution of the imaginary part DSEs. Before making a comparison of physical 
quantities we present the numerical results for the Green’s functions. In Fig.l the so 
called dynamical mass 


M{p 2 ) — G 1 {p 2 ) — p 2 (31) 

of <f> 3 theory boson is presented for various coupling strengths in both renormalization 
schemes. The infrared details are displayed in Fig.2. The dynamical mass is not 
directly physically observable since it is scheme dependent from the definition, the 
exception is the pole mass which is scheme independent and renormgroup invariant 
as well. It is interesting that there are time-like values of square of momenta where 
the propagators behave almost like free ones no matter how the coupling constant is 
strong. This happens somewhere around the point pj = 6 m 2 for OSR scheme and 
approximately at p 2 = 20m 2 for MMS scheme, which implies the physical irrelevance 
of such a behavior (Of course, there are always differences within the absorptive 
parts up which are ordinarily coupling constant dependent at these points). The 
appropriate relevance of propagator dressing is best seen when the dressed propagator 
is compared with the free one G = (p 2 — m 2 ) -1 . From the Fig.3 and Fig.4. we 
can see that the propagator function is the most sensitive with respect to the self¬ 
energy correction for threshold momenta where these correction are enhanced about 
one magnitude, while they are largely suppressed for the above mentioned values 
of momenta p 2 . Note that nothing from these things can be read from the purely 
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Euclidean approach. The results presented up to now have been calculated in the bare 
vertex approximation, the solution with vertex correction included will be discussed 
bellow. The appropriate bar vertex approximation critical coupling value is X^ r u — 
3.5 for OSR scheme and A — 5 for MMS scheme. Their different values are not a 
discrepancy but the necessary consequence of the renormalization scheme dependence. 

Furthermore, in order to see the effect of self-consistency of DSE treatment we 
compare the DSE result with the perturbation theory in OSR scheme. From the 
Fig.5 we can see that the perturbation theory is perfectly suited method when applied 
somewhere bellow the critical value of the coupling. Therefore, the main goal of our 
solution is the information about the domain of validity of given model. 

The issue of vertex improvement by the one loop skeleton diagram and its ap¬ 
propriate effect on the DSEs solution and scattering matrix is discussed in the text 
bellow. First let us note that the critical values of the couplings decrease and we 
have \dressedvertex ^ ybarevertexj 2 w j 1 j c ] 1 j s roughly valid for both the renormalization 
scheme employed. We return to the question of meaning X crit when we will discuss 
the WCM. 

To make our comparison of proposed methods more meaningful, we do not com¬ 
pare the Green’s functions but rather wee look on the scattering amplitudes M cal¬ 
culated in both renormalization schemes obtained in both truncations of DSEs. In 
Fig.6 we compare the imaginary parts of scattering amplitudes M at a given kine¬ 
matic channel. The comparison is made in the way proposed and described in the 
previous section. Henceforth, what are actually compared in this Figure are the 
Lehmann weights <r’s of the MMS scheme calculated for certain A mms an d the rescaled 
Lehmann weights ro 2 calculated for the OSR scheme with the appropriate coupling 
strength A osr — v Xmms ■ If is apparent that the lines for ‘XsmMosR^XMMSit) and 
^smM MM s(XMMS,t) f° r solutions with dressed vertices are much close each other 
then the solution with bare vertices. This statement is valid for all t for a given 
theory characterized by its coupling constant (with A mms fixed). This is true for 
all couplings X's, the only-but not so striking- exception is certain infrared excess 
for the value of couplings closed to the critical one. Of course, the worse numerical 
accuracy play the role in strong coupling. Nevertheless, we can see that when the ap¬ 
proximation improves then there is apparent signal for achieving the renormalization 
scheme independence for all the values of the coupling constant. 

In order to see aforementioned quantitative improvement we have calculated the 
appropriate deviations E n ,N = 0,1. The results for some larger value of the cou¬ 
plings are presented in the Tab.l. The corresponding difference becomes negligible 
when A decrease and approaches its ’perturbative’ value. For better orientation the 
infrared details for three choices of the coupling constants are also displayed in Fig. 7. 


3 DSEs for the WCM 

The massive WCM is given by the following Lagrangian 

£ = £ - y \ $3 + c.p. 


(32) 
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where C.P. means the appropriate counter-term part. Here we choose the second 
renormalization scheme employed in the previous section, i.e. the propagators of all 
three particles have the unit residua. All the definitions of counter-terms SZ t , Srrii, Sgi 
correspond with the OSR defined previously but now for each particle separately. 
Furthermore, we adjust the couplings to be 

to = ^rto 0 , i = 1,2 (33) 

z&l 


such that 


.913 — #23 • (34) 

The equal mass case rrii = m 2 was already solved [24] for purpose of studying 
the self-energy effect on the bound state spectrum. Here we solved the unequal mass 
case 


— = 4; m 3 = m 2 (35) 

m 2 

and compare the result with the Euclidean version of solution [16]. We restrict 
ourselves to the bare vertex approximation which is sufficient for comparison with 
[16]. Since all the derivation is rather straightforward we simply review the results. 
The renormalized DSEs in bare vertex approximation read 


GrI ip) 

= p 2 - 

m 2 i - U i{2] (p 2 ) i — 1,2 


GrI {p) 

= p 2 - 

ml-n m (p 2 ) 


n i(p 2 ) 

= *2 g 2 

J ( 2n yi ip q ) Gi iq ) % = 1 , 2 


n 3 (p 2 ) 

= 

f ( ,4 5 ] Gi ip q) Gi iq) 

(2vr) i=12 

(36) 


where the bracketed index denotes the renormalization scheme employed,and the sec¬ 
ond index labels the particle associated with the appropriate held in the Lagrangian 
(32). All the propagators satisfy the Lehmann representation with unit residuum 
and all the proper function obeys the double subtracted dispersion relation (17). 
Henceforth, the appropriate spectral weights are related through the relations 


= 


Pi(v) 


(U 


m 2 ) 2 


+ 


(u 


m7 


I (mj+m k ) 2 


da- 


(k-’) Pi 


a) 


1 2 


+ <Ji(a)pi( u) 


LU 


a 


(37) 

where for the indices i — 1,2 we have j = 2, 3; k = 3, 3 and for particle 3 the index 
j — k — 2 since it label the lighter particle with the mass m 2 . The expression for the 
absorptive parts read 
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= 


2 g 2 
(4tt) 


B(rnf, m 3 ; u) + r do; (l>(ck, ml',oj)(7z{a) + f?(a, m 3 ; u;)cri(a)) 


dad/3B(a , /?; u;)cr 3 (a;)crj(/3) 


i = 1,2 




i=lg ( 4?r ) 2 
roo 

/ dad/3B(a, /3]Lu)<Ji(a)a, 


4777,? coo 

1-- + 2 / daB(a,m%\oj)ai(a) 

to Jo 


(38) 


where we freely integrate over the whole range of positive real axis leaving the infor¬ 
mation about the appropriate thresholds and subthresholds absorbed in the definition 
of the function B. 

The above set of equations has been actually solved numerically. The main result 
for us is the appearance of the critical coupling strength A c = g^/(Anml) = 0.12 which 
rather accurately corresponds with the point where the renormalization constant 
Z 2 turns out to be negative. The appropriate dependence of the renormalization 
constants Zi is presented in Fig.8 for all three particle. The obtained critical value 
is in reasonable agreement with the one obtained by the Euclidean solution of DSEs 
system [16], where A c = 0.086, as well as with the critical value A c = 0.063 which was 
found using a variational approach [18],[19]. 

Furthermore, the existence of the critical coupling of OSR scheme can be seen 
from the analytical formula for the inverse of propagator 


G i 1 = p z - m z - 


da 


Pi( a )(p 2 


mi) 2 


0 (a — m 2 ) 2 (p 2 — a + ie) 


(39) 


which implies that for the strong enough coupling the Landau pole should appear, 
which must arise when the factor L 



(40) 


is negative (when it is just zero then the Landau pole is situated in space-like in¬ 
finite, and for the positive L this singularity never appears due to finiteness of the 
appropriate integral in (40)). For negative L the propagator cannot satisfy Lehmann 
representation at all and at least the Minkowskian treatment used in this work must 
fail. Comparing equation (40) with the definition of renormalization constant Z we 
clearly have the identification L = Z. As we have mentioned, the numerical solution 
start to fail when the condition Z — 0 is fulfilled. This statement is justified with 10 
% numerical accuracy. (We have no similar guidance for MMS scheme but we expect 
the similar appearance of the critical coupling A mms f° r this scheme as occured in 
<f> 3 theory, but the reason for the numerical failure in this case is not found in this 
paper.) 


12 



4 Conclusion 


We have obtained numerical solutions of the DSEs in Minkowski space for <f> 3 theory 
and the WCM. This suggests that the expansion of the theory around the metastable 
vacuum leads to the predicative result. Our technique allows us to extract prop¬ 
agator spectral function p(s) with reasonably high numerical accuracy. Since the 
renormalization procedure is performed analytically, it has no effect on the precision 
of solution. When the coupling does not exceed a certain critical value, then the do¬ 
main of analyticity of the propagator is the all real axis of p 2 . An attempt to clarified 
the meaning of critical coupling value was made. This suggests that it corresponds 
with appearance of unphysical singularity in the on-shell renormalized propagator. 
Consequently, the field renormalization constant (in on shell scheme) turns to be 
negative for A > A crit . 

A Dispersion relations for self-energies in bare 
vertex approximation 

In this Appendix we derive DRs for self-energies in both renormalization schemes for 
the bare vertex. The calculation is very straightforward, and in fact it represents 
nothing else but evaluation of the one loop scalar Feynman diagram with different 
masses in internal lines. 

Substituting the Lehmann representation for MMS propagators (12) the unrenor¬ 
malized II can be split like 


n (p 2 ) 
n ( b ,b)(p 2 ) 

n (M)(p 2 ) 

n (s,,)(p 2 ) 


n (6,6)(p 2 ) + 2n (M (p 2 ) + Ti( S ,s)(p 2 

'■ 18 r 2 g 2 


dq 


((p + q) 2 — m 2 + ie){q 2 — m 2 + ie) 
18 rg 2 a(a) 


dq j da 


dq / dad/3 


(q 2 — a + ie){{p + q) 2 — m 2 + ie) 
18 g 2 <j(a)a(/3) 


((p + q) 2 - a + ie)(q 2 — (3 + ie) ’ 


(41) 


where we have used shorthand notation for the measure id 4 q/( 27t) 4 = dq. Making 
the subtraction, we immediately arrive for the pure perturbative contribution (up to 
the presence of the square of residuum): 


n 


l (b,b) 


du- 


pi(b,b)(p 2 ~ m 2 ) 


Pi(b,b)(u) — 


Jim 2 (p 2 — uj + ie) (c o — m 2 ) 
18 r 2 g 2 I 


Am 2 


(4tt) 


cu 


(42) 


The most general integral to be solved is similar to the above case but with the 
physical masses replaced by the spectral variables. 
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I(p2> / d9 (0> + 9) 2 -«)(9 2 -ffl <43) 

which after the subtraction (44) and integration over the Feynman parameter x leads 
to the appropriate single-subtracted DR (44) 


hs(p 2 ) = /(p 2 )-/(m 2 ) = 


di] 


(An) 2 Jo 


dx 


poo 

n 2 


du 


(p 2 — m 2 ) 


m Z-% x) ~ <r n 2 )(p 2 — u + ie) 


du- 


P 


m 2 B(a,/3]u) 


(44) 


Jo (p 2 — u + ie) (An) 2 (u — m 2 ) ’ 
where the function B(u,v,u) is defined through the Khallcn triangle function A like 


B(u, v, u) 
A (u, c l>, v) 


\ 1/2 (u,u,v)^ / . i 

-——-—-0 (c<J - oa +/3 2 ) 

cu v 7 

(w — uj — v) 2 — 4c ov 

u 2 + u 2 + v 2 — 2uv — 2 uu — 2 uv. 


(45) 


It can be easy checked that B(m 2 , m 2 ; u) = y^l — ^-Q(u — 4m 2 ) which was already 
introduced in (42). 

The OSR scheme requires additional subtraction which is finite and henceforth 
can proceed by making a simple algebra 


h S (P 2 


= hs(p 2 ) - 

1 r x 


_d_ 

dp 2 lp 

du 


2 =m 2 ^ls (P ) 


(p 2 


m 


2\2 


(47r) 2 Jo (cu — m 2 ) 2 (p 2 — ceJ + ie) 


R(a, /?; cu) 


(46) 


To summarize the results we see that MMS self-energy satisfies one subtracted 
DR with the absorptive part npi given like 


PiriM 

Pl(fe,6) M 

Pl(b,s) (^) 

Pl(s,s) ( < ^) 


Pl(b,b) (<^) + 2pi(6 )S )(cj) + Pl(s, s ) (c<j) 

18r 2 g 2 / 4m 2 " 

(47t) 2 V a; 



da 


18rg 2 

(47t) 2 


B(a, m 2 ; 


cj)(r(a) 



R(a, /3; cj)(r(a)(T(/3) 


(47) 


while the self-energy in OSR scheme satisfies double subtracted DR with the absorp¬ 
tive part iip 2 
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Pit 2 (^0 

— P2{b,b)ipj) + 2p 2 (b,s)(<^) + P2(s,s)(pj) 


18 g 2 I Am 2 

P2{b,b)\U) 

(47t) 2 V w 

P2(b,s){u) 

= 1 da . A . 0 B(a,m 2 ;u)a 2 (a) 

J9m 2 (All)- 

P2(s,s)(v) 

f°° 18 g 2 

= / dadf3— — -B(a, /3;ou)a 2 (a)c 

J 16m 2 (47T) Z 


( 48 ) 


B Two-loop skeleton self-energy DR 

The finite two-loop integral appears after the substitution of the vertex (11) to the 
self-energy formula 


/ 


dqdk 


[( k 2 - CD + ie)((p + k ) 2 - a 2 + ie)((k - q ) 2 - a 3 + ie) 
(i q 2 - «4 + ie)((p + <?) 2 — Q(5 + ie)l 


(49) 


where all the irrelevant pre-factors are omitted for purpose of the brevity. They 
will be correctly added at the end of calculation for both renormalization schemes 
separately. The contribution is ultraviolet finite therefore we first calculate the un¬ 
renormalized result. Firstly, we parameterize the off-shell vertex by matching the first 
three denominators and consequently we integrate over the triangle loop momentum 

k 


( k 2 - + ie)((p + k ) 2 - a 2 + ie)((k - q) 2 - a 3 + ie) 


-l 


k J dxdy2y yt 2 xy + (p + k) 2 ( 1 — x)y + (k — q) 2 ( 1 — y) 
OL\xy - a 2 (l - x)y - a 3 (l - y) + ie] -3 = 


(50) 


dxdyy r 


P (1 - a) 2 /(l - (1 - a;)j/) + 0(1 - y) + 2p.g(l - a)y(l - y) 


Jo (An) 2 L j 
-cxixr/ - a 2 (l - x)y - a 3 (l - y) + ie] 


-i 


Next, we substitute x —> 1 — x and after some algebra we obtain for equation (50) 


dxdy 


Jo (47t) 2 (1 - y) 


2,o , 2 ^( 1 xy ) ^ , . 

q + 2 p.ga: + p — ---Oi _ 3 + ie 


-l 


(51) 


where we have used short notation Cfi _ 3 = + a 3 ^- We continue 

by matching equation (51) with two spare denominators in (49) by using Feynman 
variables z and u for denominators with ctq and 0 ) 5 , respectively. Then we can write 
for (49) 
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dq 


rl dxdydzdu2u p 2 


P 


Jo 47 r) 2 (l — y) L 
2 (1 — xy)xzu 


q 2 + 2p.qxzu + 2p.q(l — u ) 


1 -3 


( 1 - 2 /) 


+ p (1 — u) — Oi _5 + ie 


( 52 ) 


where we have used shorthand notation O 1-5 = Oi- 2 ,zu + 0 : 4(1 — z)u + < 25(1 — u). 
Shifting q + p(xzu + 1 — u) —> q and integrating over new q it yields: 


dxdydzduu 


1 0 ( 47 r) 4 (l - y)F(x,y,z) 




(1 — xy)xzu , .. o 

F(x, y, z) = 1 - u + ^ _ -(a^w + (1 - u)) 2 

Next we substitute u —> u where u = 0i_ 5 / F(x,y, z). Using the notation 


(53) 


U) 


CLi 


a 2 

h 


VLCL\ + ( 2-2 

m 2 6i + ub 2 

(c*i(l - x)y + a 2 xy + a 3 (l - y))z + (o 4 ( 1 - 2 ) - « 5 )j/(l - j/) 

“52/(1 -2/) 

-(1 -xz) 2 y( 1 - y) 

(1 - 2xz)y(l - y) + (1 - zy)zyz , 


(54) 


(55) 


we can write down the appropriate DR for equation (49) 


D(u;;q:i, .. a 5 ) 


U 

A 


du 


dxdydz 0 (u; — 


fSff) ep) 


/o p 2 -w + ieio ( 47 t ) 4 (1 — y) [p - W (1 - xz) 
-B + /D 


(56) 


2 A 


D = B 2 - 4AC 


ubi; B = u>b 2 — oq; C = —02 . 


Note here that spectral function (everything after the first fraction in (56) is 
always positive for allowed values of a's and it is regular function of its argument u. 
The various subthresholds are then given by the values of Lehmann variables ck’s in 
accordance with the step function presented, noting that the perturbative threshold 
is given again by 4m 2 and in that case case the result partially simplified. For 
completeness we reviewed the associated simplifications, namely: a 1 = m 2 z(l — y(l — 
2 /)); a 2 = m 2 y(l — y). Making one subtraction for the MMS and two subtraction for 
OSR scheme we can recognize that the appropriate skeleton DR for master diagram 
has the absorptive part 



(6 g) 4 
2 


H J daid(ai)Q.(u>, oq, ..a 2 ) 
1=1 ' 


(57) 
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for MMS scheme and 


p [ 2 ] (uj) = n / daia 2 (ai)Q(uj,ai, ..a 2 ) (58) 

i =1 J 

for OSR scheme, respectively. In fact it gives rise 28 various contributions to p^ (only 
12 are actually topologically independent, distinguished by the number of continuous 
Lehmann weights with the appropriate position of spectral variable in fh) All of 
them have been found numerically for the purpose of DSEs solution. 
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> 

m 

o 

II 

fe; 

bq 

En= o DV 

E n= \ bv 

Ejy = 1 DV 

A = 0.1 

0.020 

0.0039 

0.022 

0.0078 

A = 0.5 

0.076 

0.025 

0.071 

0.025 

A = 1.0 

0.15 

0.040 

0.13 

0.08 


TABLE 1. Normalized and weighted integral deviations E^ between scattering 
amplitudes calculated in the OSR and MMS scheme. Exact scheme independence 
corresponds to the case En = 0. The parameters N = 1 makes the quantity E more 
sensitive to the systematic error in the infrared domain. The denotation BV(DV) 
means that the appropriate propagator was calculated with bare (dressed) vertex. 
The function E is displayed for three cases of results presented in Figures 6 and 7. 


19 



















P [m] 

Figure 1: Dynamical mass of scalar particle in <f> 3 theory calculated in bare vertex 
approximation in the both renormalization schemes. The lines are labeled by the 
value of Xmms f° r MMS scheme and Xqsr f° r OSR renormalization scheme. 



p [m ] 


Figure 2: Infrared (threshold) details of the Fig.l. 
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Figure 3: The propagators deviations from free theory. The propagator is calculated 
in minimal momentum renormalization scheme for various A mms 



P [m ] 


Figure 4: The propagator deviation from free theory. The propagator is calculated 
in on mass-shell renormalization scheme with unit residuum for various Xqsr- 
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Figure 5: Comparison of DSE results in bare vertex approximation with the per¬ 
turbation theory result. DSE and bubble summation is compared in OSR scheme. 
Each two close lines off different types correspond to the same value of coupling 
A osr — {0.25; 1.0; 1.5; 2.2} . The lowest dashed line with A = 10 has not its DSE 
partner solution (since A > A c ). 
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Figure 6: Imaginary parts of scattering matrix calculated with propagator which 
have been obtained in MMS and OSR scheme with (dv) and without (bv) improved 
vertex. Each set of lines corresponding to the same model is labeled by the coupling 
strength A M ms- 



Figure 7: The low frequency details of the Fig.6. 
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Figure 8: The dependence of field strength renormalization constants on the coupling 
strength of Wick-Cutkosky model. The index 1-3 labels the particle. 
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